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Abstract 

A simple model of the Majorana neutrino with the see-saw mechanism is 
^ \ studied, assuming that two light neutrinos are the hot dark matters with equal 

mass of 2.4eV in the cold plus hot dark matter model of cosmology. We find 
that the heavy neutrino, which is the see-saw partner with the remaining one 
light neutrino, can be the cold dark matter, if the light neutrino is exactly 
massless. This cold dark matter neutrino is allowed to have the mass of the 
wide range from 5.9 x lO^eV to 2.2 x lO^eV. 
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Recently Primack et al. |T[ pointed out that the cold and hot dark matter model agrees 
very well with the observations of the matter distribution in the universe with the total 
density parameter = 1 and the Hubble constant h = Hq/IOO km s~^ Mpc"^ = 0.5. They 
assumed that two massive neutrinos which have nearly degenerate masses 2.4eV play a role 
of the hot dark matters. The hot dark matter, the cold dark matter, and the Baryon occupy 
20%, 72.5%, and 7.5% of the total density parameter, respectively. 

On the other hand, there are some current status for the masses and flavor mixings of 
neutrinos. The solar neutrino deficit and the atmospheric neutrino anomaly |^ seem to 
give the indirect evidences of the non- vanishing masses and flavor mixings of the neutrinos 
in the view of neutrino oscillation. In addition, recent LSND experiment seems to have 
brought the first direct evidence for neutrino masses and flavor mixings in the ^ 
oscillation. If some species of neutrinos have masses of order eV, they can be appropriate 
for the hot components of the dark matter in the cold and hot dark matter model. 

In the standard model of the elementary particle physics, three species of neutrinos are 
exactly massless and there is no particle which can be the cold dark matter. Some extension 
is needed to include the mass of the neutrinos and the cold dark matter. In this letter, we 
study the model first introduced by Chikashige, Mohapatra and Peccei [0. This model is a 
very simple extension of the standard model, which includes massive Majorana neutrinos. 

We introduce three species of the right-handed neutrinos and an electroweak-singlet 
scalar $ as new particles to the standard model. In order to make neutrinos massive, two 
kinds of the Yukawa couplings are considered. The Yukawa interaction is described by 

-^Yukawa = - 9 y ij'^i^'^ Rj " 9 Aiij'^Ri^'^Rj + h.C, (l) 

where (p is the electric-charge neutral component of the Higgs field in the standard model, 
and i and j denote flavors = 1,2,3). The Dirac and the Majorana mass terms appear 
by the non- zero vacuum expectation values of these scalar fields ((0) 7^ and ($) 7^ 0). The 
mass matrix is given by 
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rriD 



+ h.c, 



(2) 



where m/5 ( M ) is the 3x3 Dirac (Majorana) mass matrix defined by m^jj = gYij{4>) {Mij = 
gj^jij{^)), and ul and i//? are the column vectors of the three fiavors. Since the symmetry of 
the lepton number is spontaneously broken by ($) 7^ 0, a massless Nambu-Goldstone boson 
called majoron appears. Generally mo { M ) is the 3x3 complex (symmetric) matrix and 
we should diagonalize the whole 6x6 mass matrix in eq.(§). 

As the first approximation, we assume that the off-diagonal elements of the mass matrices 
rriD and M are very smaller than their diagonal elements, namely, the effect of the flavor 
mixing can be neglected. By this assumption, we can take m^) = diag [1711,7712,171^] and 
M = diag [Mi, M2, M3]. Furthermore, we assume the hierarchy between mo and M, namely 
\mi/Mi\ ^ 1. Then, the see-saw mechanism separately works on each generation. 

The six mass eigenstates are described by the weak eigenstates, ul and ur, as 

(3) 

where e = diag[mi/Mi,m2/M2,m3/M3], and z/£ = (z/^ z/^ u^)'^ and = (z/^ 1^2 ^3)^ 

light and heavy Majorana flelds, respectively. The masses for z/f and are ml ~ —{miY/Mi 

and m^ ~ Mi, respectively. 

The couplings of the light and heavy neutrinos with the weak bosons are given by 

e 
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7m(1 - Ibji^t + ^7/^(1 - 75)e'^hJ + h.c. , (4) 
17^7^(1 - 75)z/£ + {TT'af.il - l-,)evh + h.c.} + 17^7^(1 - l->)e^Vh\ , 

where £ denotes the column vector of three charged leptons, i = {e fi t)'^. Here we neglect 
the flavor mixing also in the charged leptons. The couplings of the heavy neutrinos with 
the weak bosons are suppressed by the small factor e. These flelds ui and z//j also couple 
with the majoron (imaginary part of the fleld $) Q According to eq.(|l|), these couplings are 



^ These fields and also couple with the majoron-Higgs (real part of the field and the usual 
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given by 

i^x^ - --^ X iibi^i - h5i^h + h.c.} + i7Kg°h5i^h] , (5) 

where the field x is the majoron field defined by x/ = Ini$, and g'^^ is a diagonal matrix, 
(yf^^ = diag[gi , g2 , g^] = M/ ($). In contrast with eq.(|D, the couplings of the light neutrinos 
with the majoron are suppressed by the factor e. 

Some parameters are fixed according to the cosmological model considered by Primack 
et al. In their model, the mass spectrum of the light neutrinos is constrained as m\ <^ 
m\ — iTil ~ 2.4eV. Then, the five parameters are left free in our model, gi' , g2 , gf ■, 
and m\ (<^ 2.4eV). Since it is very complicated to analyze leaving all of these parameters 
free, we consider the simplest case in this letter. We set gi = g2 = g" = gM-i therefore 
m\ = 1712 — ^3 = ^h- Now we have only three free parameters, g^^^, rrih, and The value 
of m{ is fixed as the following discussion. 

There is a fundamental constraint that the particle has to be stable to be the cold dark 
matter, namely, its life time must be longer than the age of the universe {tu ~ 10^'' s). The 
heavy neutrinos decay into the light neutrinos and the majoron through the coupling in 
eq.(^). The life time of the i-th heavy neutrino is given by 

327r 

Unless gj^j is very tiny (g^^j < 5.3 x 10^^^), two neutrinos, z/g and z/g, are unstable because 
of — ml ~ 2.4eV. Since such extremely small coupling constant is unnatural, these 
two neutrinos cannot be naturally the cold dark matters. Then, the heavy neutrino is 
the only remaining candidate for the cold dark matter. If we adopt natural value (not so 
small value) for the Yukawa coupling g^^j, the mass of the corresponding light neutrino z/f 
must be extremely small, for instance, m{ < 6.6 x 10~^''eV for g^^^ ~ 10"^. Since we cannot 



Higgs in the standard model. However, we do not consider the effects of these fields assummg that 
masses of these couplings assuming that these scalar particles are very heavy. 
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believe that such small value of the mass is explained by some mechanisms (Dirac mass mi 
must be extremely smaller than the weak scale), we assume that is exactly massless by 
virtue of some symmetries. Note that the vanishing m{, or the Dirac mass mi = 0, is stable 
against the radiative correction, since the U{1) symmetry of the phase rotation of the fields 
I'll and e forbids the generation of the Dirac mass mi. Furthermore, it should be stressed 
that only the mass difference is important for the neutrino oscillation phenomena, but not 
the absolute value of the mass. There is nothing wrong with the massless z/f. According 
to this assumption, the heavy neutrino Ui can be stable cold dark matter. In the following 
discussion, our aim is to investigate the cosmologically allowed region of the mass mh and 
the coupling constant Qj^j. 

Let us consider the dynamical evolution of the early universe in this model. At the 
very high temperature, we assume that all the particles are in thermal equilibrium. As 
the temperature cools down, some particles decouple from the thermal equilibrium at each 
specific temperature called decoupling temperature. It is convenient for considering the 
evolution of the universe to separate the matter contents into two parts. One is the 'heavy 
neutrino-majoron system' which includes three heavy neutrinos and the majoron. The 
other is the 'electroweak system' in which all the other particles are included. These two 
systems weakly interact with each other through the couplings suppressed by the see-saw 
factor e in eqs.(^ and @. We assume that at the temperature the 'heavy neutrino- 
majoron system' decouples from the 'electroweak system', and that the particles in the 'heavy 
neutrino-majoron system' are still in thermal equilibrium after this decoupling. Therefore, 
the another decoupling temperature is lower than . Here is the temperature at 
which the heavy neutrinos and the majoron no longer interact with each other. 

Note that the temperatures of these two systems are different after the decoupling at 
rpEW ^ The temperature of the heavy neutrino as the cold dark matter in the present universe 
TcDM is different from the temperature of photon at present, Tr = 2.7K. The reheating fac- 
tor aji is estimated by considering the reheating of photon caused by the charged particle 
and the anti-particle annihilation. In addition to this usual factor an, the cooling factor R 
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should be introduced in our model in order to include cooling effect by the decaying unstable 
neutrinos. As soon as the 'heavy neutrino- majoron system' decouple from the 'electroweak 
system' at , two unstable neutrinos 1^2 and decay into the light neutrinos and the 
majoron. Then, the 'heavy neutrino-majoron system' cools down and the 'electroweak sys- 
tem' is heated up, because the energy of the 'heavy neutrino-majoron system' flows to the 
'electroweak system' by the emission of the light neutrinos. Since the degrees of freedom of 
the 'electroweak system' is far larger than that of the 'heavy neutrino-majoron system', we 
can ignore the heating effect of the 'electroweak system'. The coohng factor R is defined by 



where T^w temperature of the 'electroweak system' when the heavy neutrinos decouple 
from the majoron. Since < as assumed above, the factor R is less than unity. The 
temperature Tcdm is written as 



by using the factors and R. 

Now we can write down the condition for the heavy neutrino to be the cold dark 
matter, p^h = pcdm- Here pcoM = 1-9 x 10~^GeV/cm^ is the energy density of the cold 
dark matter in the present universe based on the model considered by Primack et al. The 
present energy density of the cold dark matter neutrino is given by 



where n{T^) is the number density of the cold dark matter at the decoupling temperature 
Tq. This number density is given by 



Tcdm — o.r RT^ . 



(8) 




(9) 




(10) 



where E — Jp^ + ml, xd — T]^/mh, and /{xd) is defined by 
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Considering that the cold dark matter should decouple in the non-relativistic regime, we 
provide the upper bound of as < 1. By using the condition, p^^h = pcdm, and 
eqs.(^-(|TO|), the mass nih is given by 

mh{xD,R) = 1.1 X 10"^ a]^^ f{xD)~^ GeV (12) 

as the function of and R. 

On the other hand, the coupling constant Qj^^ is obtained from the definition of the 
decoupling temperature T^. The decoupling temperature is defined by 

nin) {a\v\)^. = H , (13) 

where (cr|f |)yx is the average value of the annihilation cross section of a heavy neutrino times 
relative velocity, and H is the Hubble parameter. For the non-relativistic heavy neutrino, 
we obtain 

by considering the heavy neutrino annihilation process, u^u^ XX- The Hubble parameter 
H is given by 

H= (15) 



^ 90 ; Mp 

where Mp ^ 1.2 x lO^^GeV is the Planck mass, and (7* is the total degrees of freedom of the 
all particles in thermal equilibrium. Note that T is not equal to T^. Since the degrees of 
freedom of the 'electroweak system' are far larger than that of the 'heavy neutrino-majoron 
system', the expansion rate of the universe, or the Hubble parameter, is approximately 
controlled only by the 'electroweak system'. Therefore, we can set T and g^^ the temperature 
T^^r and the degrees of freedom of the 'electroweak system', respectively. From the definition 
of R in eq.(^, the Hubble parameter is rewritten as 



\90 J MpR^ ^ ' 

Substituting eqs.([TO|), (p^) and (0) into eq.(p!3|) and eliminating nih by using eq.(|l2|), g^j is 



given by 

g,,{xn. R) = 5.0 X 10-^ g^' a^'/' x~^"' fi^xnY'" R-'" • (17) 

as the function of xd and R. Since we obtain both irih and g^^^ as the functions of x^i and 
i?, one hne is drawn in the mh-g^j plane for one fixed values of -R(< 1) varying xn from zero 
to unity. The allowed region of nth and g^^^ is very large, if the value of R is absolutely free. 

Next we estimate the cooling factor R by using the 'sudden-decay' approximation for 
two unstable heavy neutrinos z/g and z/g . We approximately consider that all the unstable 
neutrinos decay and disappear at once when the age of the universe is equal to their life 
time, r {= T2 — t^) 0. In addition to this approximation, we assume that the disappeared 
z/g and z/g are quickly supplied by the majoron annihilation, and the thermal equilibrium is 
recovered. The same situation is expected to occur also at the age t = 2r, 3r and so on, until 
the temperature of the 'heavy neutrino- majoron system' cools down to the decoupling tem- 
perature T^. According to these approximations, the ratio T^/T^ can be estimated, where 

is the temperature of the 'heavy neutrino- majoron system' just after quick supplement 
and is the one just before the 'sudden-decay'. 

The energy density of the 'heavy neutrino- majoron system', phm, is described in two 
different ways. Since the heavy neutrinos and the majoron are in thermal equilibrium just 
after quick supplement, we obtain 

Phm = ^(1 + ^x3) Ti . (18) 



^ The unstable neutrinos 1^2 and u-^ do not start decaying from t = 0, but t = at which they 
decouple from the 'electroweak system'. However, we can approximately set = 0, because 
r is satisfied in our final result. 



8 



On the other hand, just after 'sudden-decay' (before quick supplement), phm is given by 

P^- = + i X 3) n - X 2 X i) . (19) 

The second term denotes the loss of the energy density due to the emission of the light 
neutrinos. From these two expressions of phm-, we can obtain T^jT^ = (18/25)^/^. Since 
the 'sudden-decay' and the quick supplement recurrently occur, we obtain Thm/TEw = 
(18/25)'^/'^ at the age t = nr, where Thm and Tew are the temperature of the 'heavy neutrino- 
majoron system' and the 'electroweak system' at the same age, respectively, and n is the 
positive integer. This ratio is translated to the smooth function of the age of the universe t: 
Thm/TEw = (18/25)*/'*'^. Therefore, the cooling factor R is given by 

(-) 

where is the age of the universe at which the heavy neutrino decouple from the majoron. 
On the other hand, the definition of R in eq.(0) is rewritten as 



i?=^ = ^. (21) 



Since T^vk described by t^, by using the relation between the Hubble parameter and the 
age of the universe: 



we obtain 



By using eqs.(||), (|2T|), and (^3]), we can obtain a relation among m/j, Qj^, xd, and as 



V25y T|^(t 



Now we obtain three independent relations, eqs.(|T2D, (p!7|), and (0) for m/i, g^^, xd and t^). 
Therefore, we can draw one line in the rrih-gj^,, plane. The result of numerical calculations 
for these relations is shown in Fig.l (upper solid line for various values of X15 < 1). 



However, note that our approximations underestimate the value of R. Because the correct 
amount of the decaying i/g ^ind z/3 is clearly smaller than that estimated by 'sudden-decay' 
approximation. Furthermore we provide next decaying neutrinos by quickly supplement 
approximation, although the disappeared z/j and z/3 are not so quickly supplied. Since both 
rrih and g^^ are decreasing functions of R as can be seen in eqs.(0) and (P!7|), the upper solid 
line in Fig.l is interpreted as the upper bound of the allowed region. 

There exists the upper bound of m/^, since we assumed that all particles are in thermal 
equilibrium at very high temperature. Then there should exist a value of temperature 0) 
which can satisfy the condition as follows: 

^>l. (25) 

Here n is the number density of the heavy neutrino and is the average value of 

the annihilation cross section of the heavy neutrino times relative velocity. Considering 
the annihilation processes I'hT-'h ^i^e, U., or qq (g denotes quark) according to the weak 
interaction of eq.(^), we obtain 

^i'^\^\) = ^^[ — ] L .KIT . . <.T.i ^TfT^^TTrFT' (26) 



37r2 KnihJ Jo e^/^ + 1 (4^2 - M|) 2 + M|r| ' 
where E = Jp^ + m|, ~ 2.4eV, Gp ~ 1.17x lO^^GeV"^ is the Fermi constant, Mz is the 



mass of the Z boson, and Tz = 2.5GeV is the total decay width of the Z boson Q. The factor 
is defined by = {I3 — Q sin^ OwY + (-^3)^) where I3 and Q are the third component of the 

weak isospin (±|) and the electric charge of the final state fermion, respectively. Summing 
for all the possible final state fermions, we obtain ~ 7.3. By numerical calculation of 

eqs.(^) and (PB|), we obtain the upper bound < 22 MeV. 

Finally, we consider the constraint from the big bang nucleosynthesis (BBN). The number 

of species of the light neutrinos is constrained as N^, < 3.04 f^. The contribution of new 



^ The annihilation cross section intermediated by the majoron is far smaller than that inter- 
mediated by the gauge boson in the region of rrih and of our final result. Then, it can be 
ignored. 
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particles (three heavy neutrinos and the majoron) to the energy density in the BBN era 
(~ IMeV) have to be small enough in comparison with 

PAiv = ^ ^ AiV X (IMeV)^ , (27) 

where AA^ = max(A^y) — 3 = 0.04. The energy density of the new particles is given by 

TT^ 7 

Pnev, = ^(1 + 4X3) X {o^rRY (IMeV)^ , (28) 

1 /3 

where qTr is the ordinary reheating factor at the BBN era, qTr = [g^{iy\.eY) / g^,{T^^)^ 
Therefore, we can obtain the upper bound of R: 

R<^-^xANy a~n-' . (29) 

This provides the lower bond on the allowed region of rrih and g^^j, since both rrih and g^^ 
are the decreasing functions of R. The bound Q is shown in Fig.l as the lower solid line for 
various values of x^- 

Our final result for the mass of the cold dark matter and the coupling to the majoron is 
shown in Fig.l. The region among the dashed line, the upper solid line, the lower solid line, 
and the horizontal line of the upper bound on rrih is allowed in our analysis. The allowed 
region of rrih and gj^j covers over about five and three orders of magnitude, respectively. The 
mass matrices of the see-saw type are realized in this allowed region of ruh. The right hand 
side of the dashed line in Fig.l satisfies the constraint that the cold dark matter should 
decouple from the majoron in non-relativistic regime. Requiring the fact that the heavy 
neutrino has been in thermal equilibrium of the 'electroweak system' once, the upper bound 
on rrih (< 22MeV) is obtained. The 'sudden-decay' and 'quick supplement' approximations 
provide the upper bound of the allowed region (upper solid line in Fig.l). The lower bound 
of the allowed region is obtained by the BBN constraint (lower solid line in Fig.l). 



^ If we refer the more loose bound, N^, < 3.3 the upper bound of R becomes little larger and 
the allowed region becomes a little larger. 
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Here we would like to mention the constraint from the characteristic mass scale of the 
free streaming of the cold dark matter. The free streaming length is roughly estimated as 



where Mq is the solar mass. This scale means the lower limit of the scale of structure which 
can be formed by the effect of the cold dark matter. If we consider that the scale of the 
globular clusters (lO^M©) should be explained by the cold dark matter, more strict lower 
bound on rrih is obtained. From the condition of Mps < lO^M© and using eqs. ([T2|) and (p!7D, 
we obtain the bound nih > 6.3 x lO^eV. 

In conclusion, we studied whether the heavy Majorana neutrino can be the cold dark 
matter or not in the cold plus hot dark matter model considered by Primack et al. The 
model of the Majorana neutrino first introduced by Chikashige, Mohapatra, and Peccei was 
considered as the simple extension of the standard model. We found that if a light neutrino 
is exactly massless, the heavy neutrino, which is the see-saw partner of the massless neutrino, 
can be the cold dark matter, provided that other two light neutrinos play the role of the 
hot dark matter. Therefore, both the hot and cold dark matters are Majorana neutrinos. 
We obtained the wide allowed region in the mh-Qj^j plane by considering the cosmological 
arguments. 








(30) 



and we obtain the characteristic mass scale of the free streaming as 
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FIGURES 

FIG. 1. The allowed region of the mass of the cold dark matter and the coupling to the 
majoron g^^j. The dashed line is the line of ,t/) = 1 for various values of R. The upper solid line 
is the upper bound for the allowed region. The dots on this line correspond to = 1, 0.6, 0.4, 0.2 
from below, respectively. The lower solid line is the lower bound for the allowed region. The dots 
on this line correspond to X£, = 1, 0.6, 0.4, 0.2, 0.1 from below, respectively. The dotted horizontal 
line is the upper bound of the mass, m < 22MeV. The region among these four lines is allowed in 
our analysis. 
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